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1. Introduction 

J GLMM 

The main goal of this paper is to generalized our approach to soft interactions developed in Ref. pi], to 
nucleus-nucleus and hadron-nucleus interactions. This approach is based on two main assumptions that 

J BFKLj LI.GLR , MUQI , MV , B, K . JI fflJLK 

provide a natural bridge to the high density QCD approach (see Refs. ^-9]). i) dp = U; and n) All 

Pomeron-Pomeron interactions can be constructed from triple Pomeron vertices through Fan diagrams. 

J GLMM 

Based on the above two assumption, we have analyzed in Ref. Jl] the available data on p — p and p — p 
soft scattering so as to determine the soft Pomeron features. We obtain: 

1) Ajp = 0.35. 

2) a'jp = 0.012. This fitted value supports our input assumption. 

3) The value of the triple Pomeron vertex coupling (?3jp = 7 Ajp is small (the fitted 7 = 0.0242). 

4) Note that Pomeron-hadron (and Regge-hadron) interactions are treated in this approach phenomeno- 
logically. 
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To summarize: The data analysis [tt] confirms our input assumptions and leads to a natural matching 
between the soft Pomeron and the pQCD hard Pomeron. Indeed Ajp ~ as, 7 ~ a s and a'jp <C 1. 

In section 2 we derive the main equations governing nucleus-nucleus scattering at high energy. To this 
end we define the kinematic regions in which our formalism is applicable, 

gA 1/3 G 3 jp exp (AY) » 1 and G\ P exp (AY) 1. (1.1) [U] 

1 + Ajp denotes the intercept of the soft Pomeron, g is the vertex coupling of the Pomeron- nucleon interac- 
tion and G^jp is the vertex coupling of the triple Pomeron interaction. As we shall see, the kinematic region 
defined by Eq. ( |l~l|) is wider than the kinematic region relevant to hadron-hadron scattering. Consequently, 
we have to go beyond the traditional Glauber- Gribov eikonal approach. The result we obtain suggests that 
it is necessary to take into account the inelastic Glauber correction which can not be expressed in terms 
of the nucleon- nucleon scattering amplitudes. In the wide range of energies where a'pY <C R 2 A , the scat- 
tering amplitude for the nucleus-nucleus interaction does not depend on the details of the nucleon-nucleon 
interaction at high energy. In the formalism we present, the only (correlated) parameters that we need to 
know are ^ffj^'ip^ an d g- Ajp, G3JP and g were obtained from the data analysis of proton-proton soft 
scattering Ji]. Our basic dynamical assumption is that a'jp = 0, which is supported by the fitted value of 
a'jp = 0.012 obtained in Ref. [lip Since the fitted values of G3JP and a'jp are small, Eq. (|l7l|) is valid over a 
wide range of energies, including the LHC energy. For the sake of completeness, we discuss in section 3 the 

main equations for hadron-nucleus interactions that have been derived in the kinematic region of Eq. (11711) 
J SCHW.'BGLM m — v 
in Refs. Jl 1^ 12] . In section 4 we adjust the general formulae of sections 3 and 4 to the specific approach 

of Ref. jlp "Section 5 is devoted to a comparison of our results with the experimental data, mostly on the 

nuclear modification factor in the RHIC range of energies. Predictions for LHC energies are presented and 

discussed. In the conclusions we reflect on the physical meaning of our approach and its relation to the 

jMy 

Color Glass Condensate (CGC) model 



2. Equations for nucleus- nucleus collisions 

JGRIBRT J COL.SQFT.LEREG 

In the framework of the Pomeron Calculus piu] (see also Refs. [113-lh]) there are two different kinematic 
domains whereone can develop a theoretical approach for nucleus-nucleus scattering. 

In the first domain we consider 
g l92 J d 2 b'd 2 b» S Al (b")S A2 (b-b')P(Y,b"-b') = g l92 j d 2 6"P(Y,6") J d 2 b' S Al (b') S A2 (b-b') (2.1) [krT 



oc 



gxg^A^A^ 3 (R Al +R%)e A - Y « 1; 



9iS Ai (b)G 3]P e^ Y oc giGsjpA^e^ « 1; G 2 31P e^ Y « 1. 



set 



In this kinematic region the main contribution stems from the diagrams of Fig. [11 Summing these diagrams, 

] GLAUB.GRrlA 

we obtain the Glauber- Gribov eikonal expressions describing nucleus-nucleus scattering |16, 1/]. Tn this 
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A, 



glset 



Figure 1: The full set of the diagrams for Glauber-Gribov approach which contribute to the scattering amplitude 
in the kinematic region of Eq. 



approach the nucleus-nucleus amplitude can be derived from the knowledge of the nucleon-nucleon cross 
section. Specifically, we get 

N d (Y;b) = i(l - exp(-^ gi g 2 J d 2 b' S M {b') S M (b -b') Jd 2 b"P(Y,b" 

= U ~ exp(~g l92 j d 2 b'S Al (b') S A2 (b-b') e AY ^y (2-2) [gg] 

In the second kinematic region 
g l92 j ' d 2 b' d 2 b" S Al (b")S A , 2 (b-b')P(Y,b" -b') = g l92 J d 2 6"P(y,6") j d 2 b' S Al (b') S M {b-b') (2.3) [KR2 

oc g 1 g 2 A\ /3 A 1 2 /3 (R 2 Al + R 2 M ) e A - Y > 1; 
giS Ai (b)G 31P e A *> Y oc 9i G^A 1 / 3 e AlpY « 1; G 2 jpe A]pY « 1, 

l aetset 

it is necessary to consider the the more complicated set of diagrams presented in Fig. |2| We conclude 
that the knowledge of the nucleon-nucleon amplitude is not sufficient to enable a calculation of nucleus- 
nucleus scattering. In addition to the above, we also need to know the structure of the Pomeron-Pomeron 

r tset 
we assume that 

only the triple Pomeron vertex contributes to multi Pomeron interactions. 

As stated, the main goal of this paper is to calculate the nucleus-nucleus amplitude in the kinematic 
region of Eq. ( ^"j ). We start with the MPSI approach 

2.1 MPSI approach 

The MPSI approach is based on the observation that for a Pomeron with a intercept larger than 1, i.e. 

fcapsiset 

Ajp > 0, the main diagrams of interest have the form of Fig. y with an arbitrary y which is of the order of 
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netset 



Figure 2: The full set of the diagrams that contribute to the scattering amplitude for the kinematic region Eq. (Eg). 



J GLMM.LMP 

y. We have discussed this approach in our previous publication Pl7T9p The method of calculating the sum 

topsiset I MPS I 

of Fig. ^ diagrams is based on t-channel unitarity adjusted [trSj to the summation of Pomeron diagrams. 



The general formula has the form 



N d (Y) - ^ (-1) -y 7 tz— p \y r =o ^— t h R =o ( 2A ) l MPSI l 



n=l 



- 1 - { exp (" 7 44) """^-^ *""(*•*)} i*-,*-- 

7^'* are related to projectile and target, respectively, and 7 is the amplitude for 'wee' partons (colorless 
dipoles in QCD) scattering at low energies. ]\f MFA denotes the amplitude, which for a nucleus has the 
form 

N% FA (y, b) = 1 - exp [-gt S Ai (b) N™ FA ) . (2.5) [jjpsil] 

Here, N^f FA denotes the sum of Fig. Jcfm grams where 51527 = 9W2- I* 1 the case of a'p -C 1, i.e. for 
energies a'jpY < R 2 A , this sum is 

N MFA _ lRe^iY-y) 

where denotes the scattering amplitude at low energy of the 'wee' parton (colorless dipole in QCD) 
with the target. 

Using a generating function for Laguerre polynomials (see Ref. jsO] formula 8.973(1)), 

(1 -*)-«-! exp (-^) = ^L%(x)z n , (2.7) [jgsis] 

^ Z lJ n=0 
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mpsiset 



Figure 3: The set of diagrams that contribute to the Figure 4: The set of fan diagrams that contribute 
scattering amplitude in the 
the kinematic region Eq. (E 



scattering amplitude in the MP SI approximatlig gjjggt I to Njf F 



312 



we obtain for Eq. (^.tj) 



N MFA (Y-y; 7R ) = - £ (&) (-^-M)" . 



(2-8 



MPSI4 



n=l 



For simplicity we have omitted the b dependence which is easy to include, as each en should be multiplied 
by S A M 



13 



Using Eq. ( |^.y[ ) and Eq. ( |2.4fj we have for the scattering amplitude, 



N 



MFA 



(2.9) 



MPSI5 



n=0 



Introducing n! = £ n e ^ d£, we can re-write Eq. ( pj.y) ) m the form 



15 



N MFA (y . 



(2.10) 



MPSI6 



n=0 



J GLMM 

70 = G3JP / Ajp is obtained from the data analysis of of soft proton-proton scattering ft] . Using formula 
JRY , 



8.976(1) of Ref. 



n! z 



n L%(x)L%(y) (xyz)-* ( 



n=0 

we derive the final result 



T(n + a + 1) 



l-z 



exp — z 



x + y 
1 - z 



In 2 



l-Z 



(2.11) 



SUML 



N™/ A (Y) 
where 



e e 1 + ^T(y) 



(2.12) 



MPSI7 



T(y) = 7e 



(2.13) 
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y 



a) 



y 




b) 



mpsieqset 



Figure 5: The diagram set that contributes to the scattering amplitude in the approximation in which only a 
merging of two Pomeron into a single Pomeron has been taken into account in the kinematic region of Eq. (fc"ljf) . 



Recalling that the function should depend on <7i<727) on qij or on q2j, one can see that Eq. ( p.l^ ) reduces 
to a simple and elegant formula for the case where we have giT{Y) ~ 1, gigiTiY) > 1 and T(Y) <C 1. 
Indeed, after integrating over £ we get 



A itk (Y; b) =1 - exp \ - 



~gfg k T(Y) 
1 + T(Y) \gt + g h ) 



(2.14) 



MPSI8 



For Eq. (\i.l4 j it is easy to write the expression that takes into account the correct impact parameter 
behavior. It has the form 



A(Y;b) 



i 1 



1 



exp 



2;7 



d Z b' 



9iS Al lb') g 2 S A2 (b-b 1 ) T(Y) 



1 + T(Y) ~ gi S Al [V) +g 2 S A2 (b-b' 



(2.15) 



MPSI91 



2.2 The complete set of equations 

We start the derivation of the complete set of equations beyond MPSI approximation assuming that the 
merging of two Pomeron into one Pomeron is equal to zero. In this case we need to sum the 'fan' diagrams 

ppsieqset ImpsjLeqset 

of Fig. |af-a which we can do using the following equation (see Fig. |t>[-b) 

Gmfa(v) = 92e^ y - G 3]P j V dy'e A ^y-y^G 2 MFA {y'), (2.16) 
which can be rewritten in a differential form 

d ° M ^ {V) = AjpG MFA (y) - G 31P Gl IFA (y'), (2.17) 



CS1 
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a) 



b) 



f ulleqset 

Figure 6: The set of the diagrams that contributes to the scattering amplitude in the approximation in which both 
the decay of a single Pomerons into two Pomerons, and^the merging of two Pomerons into a single Pomeron have 
been taken into account in the kinematic region of Eq. 




a) 



y 




b) 



mpsieqlset 



Figure 7: The set of diagrams that contribute to the Greem function G (y, y') in the approximation in which £mry 
the merging of two Pomerons into a single Pomeron has been taken into account in the kinematic region of Eq. 



with the initial condition 

G M FA(y = 0) = g 2 S A2 [b-F/) . (2.18) [csJ 

The solution to this equation is 



G M FA{y) = r 1 - t (2.19) ^S4 



1+52^ (e A ^ - 1) 

jsieqlset 





lulleqset 


mm 


and Fig 



feulleqset 

We need to sum the diagrams of Fig. B-a to obtain the exact two nuclei irreducible amplitude. We 

, Iftlset 



recall that the Glauber- Gribov formula (see Fig. 0j sums all two nuclei reducible diagrams given by the 



-7- 



exchange of a single Pomeron. The equation for the exact amplitude is illustrated in graphic form in 

j fiilleqset 
Fig. ^-b and has the form 



Gexactiv) = G M FA{Y,y) - G 3iP / dy' G 2 exact (y') G M FA(y',y)- 

Jy 



(2.20) 



CS5 



CS7 



S tapsieqlset 
The diagrams for Gm fa (y > y) are shown in Fig. Lf-a and they can be summed using the equation shown in 
sieqlset 
b. However, for the simple case of a Pomeron with ajp = we can use the property of the propagator 

for Gmfa (y) 

G M FA{y',y)G M FA{y) = G MFA (y'), (2.21) [ess 

which leads to 

G MF A(y',y) = G M FA(y')/G M FA(y). (2.22) 

Using this solution we can rewrite Eq. ( |2.2C| ) in the form 

G exac t(y) = G (^p^y) ~ ° 31P GMFA(.y) l ^ G * Xact ^ ° MFA ^ ' ^' 23 ^ 

which can be written in a differential form as 

dG exac t{y) dlnG M FA(y) 



CS8 



G exact (y) - G 31P Gi xact (y). (2.24) ^S9 



dy dy 
The corresponding initial condition is 

G exact (y = Y) = gi S Al (b') . (2.25) 



CS10 



Using the solution to Eq. (|2.24[) we can write the scattering amplitude which will differ fromEq. (fZ72|) by 
the replacement P(Y) — > G exac t [Y,b',b — b'\ . It has the form 



N el (Y;b) = i(l - exp(^-i J d 2 b' G exact (v-b' ,b - fj\\ . (2.26) 



GGE 



2.3 MPSI solution 



We have not found the general solution to Eq. (12.24), but it is easy to demonstrate that within the MPSI 

o MI 3" J [MPST8 

approximation this equation leads to the scattering amplitude of Eq. (12. 14 First, we notice that the set 

Jmpsisol ICSj 



of 'fan' diagrams in the MPSI approximation (see Fig. ||) is different from Eq. ( |2.l6| ) since, starting from 
rapidity y" the Pomerons cannot split into two Pomerons. Therefore, the number of Pomerons at y = y" 
is the same as at y = 0. Having this in mind, we obtain a solution for Gmfa (Y,y"), 

e A P (Y-y") 

Gmfa (X, y" ) = ff • (2.27) 



SDLMPSI1 



- 8 - 




One can check that this equation sums the di- 

^psisol _ jaOLHPSIl 



agrams of Fig. H by expanding Eq. Q2.2y) with re- 



spect to f^L92 e A P Y 



On the other hand, we 
can use Eq. « substituting N MFA (Y - y» j7 £) 
given by Eq. grand taking N MFA {y\ 1 t R ) = 1- 
exp (— e Ajpy ) . The second observation is that Eq. 
degenerates to 



dG, 



exact 



(v) 



dy 



A F G eX act(y) - G 3] p G 2 exact (y) 



with the initial condition 



G, 



exact 



(y = Y) = Gmfa (Y, y" ) 



(2.28) 



(2.29) 



mpsisol 



Figure 8: The set of the diagrams that contribute to 
Gmfa in the MPSt approximation. 



One can see that such a solution has the form 



G 



exact 



(Y) 



1 + §F (9i + 92) 



(2.30) 



which coincides with Eq. 



Substituting Eq. (fc.3U| ) into Eq 



since gi g 2 e 



AY 



9i 92 T (Y) and (G 3 jp/ &p(gi + 92) 



h + fc)T(Y). 



we obtain that Eq. ( gTT5p presents the scattering amplitude. 



SDLMPSI2 



SDLMPSI3 



SDLMPSI4 



3. Equations for hadron- nucleus collisions 



h 




A 





-f* — 






— * — 






— > — 




+ 


— > 






hAsct 





Figure 9: The set of diagraimthat contribute to the scattering amplitude of hadron- nucleus scattering in the 
kinematic region given by Eq. ( 3TT }T 



For hadron-nucleus collisions we have only one kinematic region (Eq. 



in which we calculate the 

scattering amplitude. This region is similar to the second nucleus-nucleus scattering kinematic region (see 
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giS Al {b)Gupe^Y oc g . G 31P A] /3 e ApY « 1; 

Gl P e AFY « 1. (3.1) [KRHA 

In this kinematic region, the hadron-nucleus scattering amplitude can be written in an eikonal form in 

i srjHW Mset 
which the opacity $7 is given by sum of the 'fan' diagrams [111] (see Fig. 



A hA iY.b) - Ml - ox P (- ^ A f ;6) ')). (3.2) ^ 



with 

qhgG enh (y) S A (b) 

ft hA (Y;6) = y^. (3.3) [M2] 

1 + gG enh (y) S A [bj 

Using Eq. ( pfJ[) and Eq. dHlaj) , we obtain that 



2 



= J d 2 b(l - exp{-n hA (Y;b))). (3.4) [HA3] 

J BGLM BDRY 

The processes of diffractive production have been discussed in Refs. pi 2, 22]. 
4. Main formulae 

In this paper Eq. ( ^.lo[ ) replaces the Glauber- Gribov eikonal formula to describe the experimental data. 
However, we need to adiust this formula to our description of hadron-hadron data given in Ref. [tTTTn this 
paper we use two ingredients that were not taken into account in Eq. ( ^.lop : 

1) A two channel Good- Walker model [El] which is exclusively responsible for low mass diffraction. 

2) Enhanced Pomeron diagrams that lead to a different Pomeron Green's function. This mechanism is the 
main contributor to high mass diffraction. 

In the two channel model we assume that the observed physical hadronic and diffractive states are 
written in the form 

if) h = a*i +/3* 2 ; V>D = -P^i + a^ 2 , (4.1) |MF1 

where a 2 + ft 2 = 1. Note that Good- Walker diffraction is presented by a single wave function tpjj. In our 
initial approach in which the Pomeron interaction with a nucleus proceeds through an elastic scattering 
with a single nucleon, we need to replace gi by a 2 g^ + (3 2 g^ 2 \ gf^ which denotes the vertex of the 
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Pomeron interaction with nucleus 1 of the states that have been described by either the wave functions 
or ^2- Since g\ = g<i we can simplify Eq. ( j^Tlj: ^"replacing g\ and gi by 9\ = 92 = 9 = ct 2 g^ + f3 2 g( 2 \ 

In the framework of our approach we consider G^jp exp (ApY) <C 1 and, therefore, we can use the 
Pomeron's Green function written in the form 



G (Y) = e A,pY . (4.2) 

I GLMM 

In Ref. |1] we sum all enhanced Pomeron diagrams. This leads to the replacement of the 'bare' Pomeron 



Green function, G (Y) = exp (AjpY), by the Green function that sums the enhanced diagrams 

1 \ 1 „ / 1 



lG{Y) 



G enh (Y) 



1 



exp 



r o, 



(4.3) 



T(Y) J T(Y) ' V"' T(Y) 

r (0, x) is the incomplete Gamma function (see 8.350 - 8.359 in Ref. [fjO])) and T (Y) is given by Eq. ( |z.l3 ). 
Finally, we have the main formulae for the total and inelastic nucleus-nucleus cross sections 





> \ 


H 


1 — exp < 

V I 



-r / d 2 b'- 



~gS Al [b')~gS A Jb-b')G enh {Y) 



G in {Ai + A 2 ;Y) = 
with ~g = a 2 ~g^ +[3 2 g( 2 \ 



d b 1 — exp 



d 2 b'- 



l + G enh (Y) \ps Al [b')+gS A2 [b 
~gS Al (b')~gS M (b-f) G enh (Y) 



l + G enh (Y)g S Al [b')+S A Jb 





Using the AGK cutting rules 

1C 



we obtain the formula for inclu- 
MJDI 

sive production (see the Mueller diagrams p4] for the process shown 



in Fig. \wfi~. The general formula for the inclusive cross section has the 
form 



1 



a in (Y) dy 



da 



a^ A {Y-y,A x )a^^{y,A : 

<Tin(Y) 



(4.6) 



Using Eq. jT.^j ) and Eq. (f^tf) one obtains 



incl 



a 



MFA 



d 2 b- 



gG enh (y)S A (b 



Figure 10: The set of the diagrams 
that contribute to the inclusive pro- 
duction of hadrom in ion-ion colli- 
sions. 



1 + gG enh (y) S A (b 



(4.7) 



MF11 



MF2 



MF3 



MF4 



MF5 



MF6 
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5. Comparison with the experimental data and predictions 



5.1 Nucleus-nucleus collisions 



The new results presented in this paper are given in Eq. Q4.4|) and Eq. ( |4.d|) . As noted they can be re- written 
in a Glauber-like form, 



1 



a^S A Ab')S A2 (b-b'\ G enh (Y) 



o tot (Ai + A 2 ; Y) = 2 / d z b 1 - exp { -- I d z b 

- I 1 + G enh {Y) [g S Al (V) + ~gS M (b - b>) 



a in (Ai + A 2 ;Y) = / d z b 1 - exp 



d 2 b'- 



ag N S A Ab')S A2 [b-b') G enh (Y) 



l + G enh (Y)~g S Al [V)+S A Jb-b> 




;(5.1) 
(5.2) 




(riot and o"^ are the total and inelastic nucleon-nucleon cross 
sections. Note that the two sets of nucleus-nucleus cross section cal- 
culations are only mildly different. A slightly larger local difference is 
observed in a restricted region of higher impact parameters. Indeed, 
for b <C R A the partial amplitude for nucleus -nu cleus scatt ering in 
either the Glauber-Gribov approach or in Eq. jrtf and Eq. 
very close to 1. At large values of b, Eq. 



37, is 

Tip and Eq. (pTzj), as well as 
the corresponding Glauber-Gribov cross sections, approach the same 
limit. A small difference between the two approaches we have studied 



St 



Figure 11: The ratio between the 
inelastic Au-Au cross sections, cal- 
culated from the exact formula, the 
Glauber-Gribov eikonal approach, 
as a function of the impact param- 
eter (b). 



tgpy> be detected only at b w 2R A (see Fig. pT±[>7 Note, in addition, 
that the contribution to the inelastic cross section in this region of b is 
very small (about 2% at RHIC energies). 

Our results suggest that sensitive observables in nucleus-nucleus 
collision are the number of participants and the nuclear modification 
factor (NMF). It is well known, that the number of participants N par t 
for Ai — A2 scattering is equal to 



N p art(b) 



A, at (b) + A 2 at (b) 



a 



A!A 2 



2 (b) 



(5.3) 



where a in % are given by Eq. (pf) and af n x (b) by Eq. (§gj. For Au-Au collisions at RHIC energies (W= 
200 GeV), the b = ratio gS§ 

N party") 



a 



NN 



'tot 



.NN 



m — v IGLMM, 
0.93 with g Gsp/Ap obtained from Ref. [ttj! In this estimate we take 

42 mb and for S A (b) we use the Wood-Saxon parameterization [Iz5], 

+00 



S A (b) 



dz 



^/z 2 + b 2 -R A ' 

1 + e h 



(5.4) 



with po = 0.1711/ fm' 

r NN ~NN _ 



NN 
'tot 



Ra = 6.39 fm and h 
36 mb. Note that, a%f f 



NN 
in 



a '">' - a el' ~ a diff = ^mo. iNote tnat, a%y f = 2a A s Y + a dd 



than the qualitative estimate of 42 mb quoted earlier (see Ref. 
We obtain N part (b = 0)/N^£ t (b) = 0.90. 

The NMF is defined as 



0.53 fm for Au. In this calculation we took a\. 
\ \ „nn _ 'p^g qJqqvq estimate is more reasonable 
) • ^mU nas been taken from Ref. Jz/f. 



Raa 



d 2 o-A 1 A 2 
dyd 2 p± 

d?(TNN 
dyd 2 pj_ 



N M ll 



d 2 N A a 2 
dyd' 2 p± 
d 2 NN N 
dyd 2 p± 



(5.5) 



NMF 



where N is the hadron multiplicity. In the last equation both d , N ^ lA2 and d , are measured experi- 

i J i dyd z p± dyd z p± r 

mentally, while the number of collisions defined as N co u = Aa^ n N /af^ A2 has to be calculated. As we have 
observed for the case of a nucleus-nucleus collision, the Glauber-Gribov approach to N co u gives the same 
result as the correct formulae, Eq. (fojr and Eq. m 4) . 



Using Eq. (g7 
Raa = 



and Eq. 



d 2 b- 



we calculate the value of Raa, 



d z b- 



S A2 (b) 



IGLMM 



1 + 9 G enh (Y/2 - y) S Al (b) J \J " " 1 + g G enh (Y/2 + y) Sa 2 (b) 1 ' 

J GLI 

In Table 1 we display the estimates of Raa f° r Au-Au collisions, calculated in the approach of Ref. [ttJT 
In this approach the contribution of the triple Pomeron interactions are relatively small, leading to a 
1 — 2mb contribution, which is just a fraction of the calculated single inclusive diffractive cross section. 
In the above, gG^jp/ Ajp = ngj (n=l,2,3) and 7 = 70 is the value obtained from the data analysis of the 



NMF1 



7/y 








1 


2 


3 


4 




inclusive 


centrality - 10% 










7o 


0.61 


0.59 


0.60 


0.57 


0.52 


0.46 


2 70 


0.42 


0.39 


0.41 


0.385 


0.345 


0.30 


3 70 


0.31 


0.28 


0.305 


0.28 


0.25 


0.218 



Table 1: Inclusive Raa for Au-Au collisions at W = 200 GeV. 
I GLMM „ 

proton-proton scattering in Ref. Y = ln(s/so) and s = W , where W denotes the center mass energy. 
Recall that Raa does not depend on p±. However, we can trust our approach only at small values of p±. 

In Table 1 we checked the single diffraction cross sections initiated by triple Pomeron interactions 

I GLMM 

corresponding to 7 values ranging from 70, obtained in Ref. [Il], to 370. Fitting the data with 7 fixed at 
either of these values we conclude that the variance of the overall x 2 /d.o.f. we have obtained is not large, 
i.e. the data is not very sensitive to the value of 7 within the above range. We believe, therefore, that the 
estimates presented in the Table 1 are instructive for obtaining an approximate value of Ra,a- For NMF 
at fixed centrality we have used Ref. p7J relations between the centrality cuts and the essential impact 
parameter region. Note that the corrections to N co u using the correct formula are small. One can see 
from Table 1 that Raa, m the centrality region (0 — 10%), is only slightly suppressed in comparison to the 
inclusive NMF (see also Eq. (p!q) ). 
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IGLMM 



kMl BRAHMSO , BRAHMS 1 , BRAHMS2 , PHENIX1 

In Fig. |l2| we plot the data for inclusive Raa (see Refs. [12S = 3I]J as compared with our predictions for 

J GL. 

W = 200 GeV. We checked two values of the triple Pomeron vertex: one which, is taken from Ref. [ttj 
(7 = 7o)> an d the second is 3 times larger. There are two different interpretations of the physical meaning 
of our results: 

1) A traditional one, in which Pomeron calculus is responsible for the structure of the initial partonic wave 

function of the fast hadron, and/or nucleus. Therefore, we need to divide the experimental values of Raa 

by the calculated NMF, and explain this ratio ^j^g^y 07 ^ by accounting for the final state interactions 

such as iet quenching, energy losses and so on [tst>, 3ot 

L ' J i glr Muqi.MV 

2) In the interpretation which we follow, the Pomeron calculus initiates the Color Glass Condensates [14-6] 

in the region of large distances. In this case, it gives the correct normalization ofRAA m the region of small 
P± *C Q s - Qs is the saturation momentum. The ratio R e f^ L /R t ^ ry can be interpreted as originating from 
two possible sources: a proper account of the transverse momentum structure of the jparton densities in 

I KHNA.KLN.KKT J ELOSS' BDMS 

the saturation region (see Refs. [kV,3z,3o]), and/or the final state interactions [136,36]. . 
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Figure 12: The experimental value of the Raa (the 

& 1 JPHENIX1 v 

data are taken from Ref. pi]) at W = 200 GeV. The 

horizontal lines correspond to two different valu A>A2 

the triple Pomeron vertex. The value of 70 as w glL. 

as other parameters such as g is taken from Ref. [t 
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Figure 13: The prediction for the inclusive Raa at 

the LHC energies. The value of 70 as well as other 

b ' I GLMM 

parameters such as g is taken from Ref. pj. 



5.2 Hadron-nucleus collisions 



The NMF for proton- nucleus collision is defined by the same expression as for nucleus- nucleus, 



d 2 (TpA d 2 N p 4 
dyd 2 p± dyd 2 p^ 
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2 


3 


4 


5 


7o 


0.9 


0.71 


0.66 


0.60 


0.53 


0.46 


0.39 


270 


0.81 


0.53 


0.495 


0.44 


0.38 


0.32 


0.26 


37o 


0.73 


0.415 


0.36 


0.30 


0.25 


0.20 


0.16 



Table 2: N 



t lie 



R P A (Y, y) N^ A eory /NpG as function of rapidity for different values of G^p. 



This ratio can be calculated from Eq. 



RpA (Y,y) 



a q~ p a 



which lead to the following equation 
gGenh(y) S A lb] 



(Pap 

dyd' 2 % 



d 2 b- 



1 + gG enh (y) S A [b 



(5.8) 



where Y = ln(s/so) and y denotes the rapidity of a produced particle. However, N co u for proton-nucleus 
scattering turns out to be different from the Glauber-Gribov approach estimates that have been used by 
the experimentalists. It means that we have to calculate 

R e Jl (Y) = RpA (Y, y) N^/Ntf, (5.9) 

where 

a pp 

A — — : — . (5.10) 



N theory (y) 



gG mh (Y)S A {b) 



For RHIC T energy range we can replace N pA eory /N pA by a, 



imation pTfT a pp (GG) = a tot - a ei 



pp (theory) 



in < 
O~tot 



gG mh (Y)S A (b) 

PF '(theory) Jo( n P 
42 mb, in Glauber-Gribov approach calculations is improved by 

IGLMM 



As we saw, a common approx- 



I GLMM ' 

o~ei — °~dif f{GW) = 36 mb, which is suitable for the approach of Ref. |1]. Recall that 

I GL 

a diff = 2o~sd + GW denotes a contribution initiated by the Good- Walker mechanism, see Ref. [Hj 

theory JJ^GG 



NMF4 



MNF4 



IGLMM 



for details and references. This difference induces a correcting coefficient N. 



pA 



36/42 = 0.86. 

However, in the framework of our approach we need to modify the total and quasi-elastic cross sections 
due to the interaction with a nucleus rather than a nucleon (see Fig. |l4[). The correct inclusion of this 
interaction leads to 



(7 in (P + A; Y) 



d b ( 1 — exp 



(5.11) 



MNF5 



a 



pp 

tot' 



S A (b) 



(l+gG enh (Y) S A (b)) 



«f + 



(7 



PP 



S A {b) 



d ^f'(l+gG enh {Y)S A {b)f 



Table 2 presents our calculations for N^ y /Ngf v As we have discussed, we can only trust these 



cou eft . ERAHMS1 

estimates at low transverse momenta. Comparing R pA (Y) with the experimental data (see Ref. [EZ9]j5 
one can see that we reproduce the low transverse momenta data well. As in the case of nucleus-nucleus 
scattering, we believe that our estimates give the correct normalization for the NMF at low p± , leaving the 
explanation of the p ± dependence of R vA to the dependence of the parton densities in the saturation region 

JKKT P I EL0gS_"BDMS 

and beyond (see Ref. [133] ) , and the jet quenching and energy loss in the final state pb, 361. If we divide 

eff J Cronin 

all data by our calculated R pA , the p± dependence of the NMF resembles the ordinary Cronin effect [134] 
which approaches unity at large p±. 
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6. Conclusions 






modn 



Figure 14: The modified nucleon-nucleon inelastic cross 
section in our approach. 



The main theoretical results of this paper are the 
formulae for the total and inelastic cross sections of 
nucleus-nucleus scattering in the kinematic region of 
Eq. ( 1^81 ) . These formulae together with equations 
for hadron-nucleus total and inelastic cross section 
and for inclusive production, enable us to make com- 
parison with the experimental data at RHIC ener- 
gies. This comparison suggests that the effect of the 
initial partonic wave function can explain the essen- 
tial part of the nucleus modification factor. Conse- 
quently, jet quenching and energy losses are respon- 
sible for a comparatively small part of the nuclear 
n | 

suppression. Since we developed our approach based 



on soft Pomeron calculus, we discuss how our ap- 
proach is related to other approaches on the market. 
There are two alternative points of view on this subject. In the traditional one, the soft Pomeron is a 
separate issue. The estimates based on soft Pomeron approach provides the information on the partonic 
wave function in the initial state. Therefore, in such an approach, the difference between the experimental 
data and our estimates should be explained by the interactions in the final state. We support the second 

I GLMH 

point of view in which the soft Pomeron approach that has been developed in this paper and in Ref. [ttjj 
is a natural generalization of the Color Glass Condensate approach, and it provides the normalization of 
the NMF at long distances. Indeed, in CGC approach the relation between N par t = cS±Q^ enters with 
the coefficient c which could only be determined from numerical simulation. We firmly believe that our 
approach suggests an alternative method to determine the numerics of CGC. The p± dependence within 
this interpretation , is correlated with the p± dependence of the partonic densities in the saturation region 
and beyond. It depends, as well, on the final state interactions, where we predict only a slight suppression. 

Using our approach we predict the NMF at LHC energies and we hope that this prediction will be 
useful for ion-ion interactions at the LHC. 
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